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^ We discuss models of interacting magnetic impurities coupled to a metallic host, which show one 

. or more boundary quantum phase transitions where the ground-state spin changes as a function of 

' the inter- impurity couplings. The simplest example is realized by two spin- 1/2 Kondo impurities 

; I ^ coupled to a single orbital of the host. We investigate the phase diagram and crossover behavior of 

PL(. this model and present Numerical Renormalization Group results together with general arguments 

' showing that the singlet-doublet quantum phase transition is either of first order or of the Kosterlitz- 

Thouless type, depending on the symmetry of the Kondo couplings. Thus we find an exponentially 
\l ' small energy scale within the Kondo regime, and a two-stage Kondo effect, in a single-channel 

situation. Connections to other models and possible applications are discussed. 
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The physics of a single magnetic impurity embedded 
in a metal - known as the Kondo effect - is a well stud- 
ied phenomenon in many-body physicsu. The low-energy 
physics is completely determined by a single energy scale, 
the Kondo temperature Tk, and the impurity spin is fully 
quenched in the low-temperature limit, T <^ Tk. 

Systems with more than one impurity are considerably 
more complicated because the interaction between impu- 
rity spins competes with Kondo screening. This inter- 
impurity interaction, which can lead to a magnetic or- 
dering transition in lattice models, arises both from di- 
rect exchange coupling /, which we explicitely include in 
the following discussion, and from the Ruderman-Kittel- 
Kasuya-Yosida (RKKY) interaction mediated by the con- 
duction electrons. 

The simplest systems containing this competition be- 
tween Kondo effect and inter-impurity interaction are 
two-impurity modelgH^Q which have been extensively 
studied in recent years. Generically, two different regimes 
are possible as function of the inter-impurity exchange /: 
for large antiferromagnetic / the impurities combine to 
a singlet, and the interaction with the conduction band 
is weak, whereas for ferromagnetic / the impurity spins 
add up and are Kondo-screened by conduction electrons 
in the low-temperature limit. 

For two spin-1/2 impurities, each coupled to one chan- 
nel of conduction electrons (total number of two screen- 
ing channels K — 2), there is perfect screening as T — > 
in both regimes and the ground-state spin will be zero. 
It has been shown that there is no quantum phase tran- 
sition as / is varied in this two-impurity Kondo model 
in the generic situation without particle-hole symmetry 
(whereas one finds an unstable non-Fernai -liquid fixed 
point in the particle-hole symmetric case)LrEl. However, 
the finite-T crossover behavior is significantly different 
in the two regimes: for large antiferromagnetic /, the 
singlet is formed at an energy scale of order /. In con- 
trast, for ferromagnetic / where both spins have to be 
quenched by the Kondo effect, the two screening channels 
are formed by the even and odd linear combiiLations of 
the conduction electrons at both impurity sitesB. These 



two channels have generically different densities of states 
leading to two different Kondo temperatures. Upon low- 
ering the temperature the two spin-1/2 degrees of free- 
dom are "frozen" at theseJjwo distinct temperature scales 
(two-stage Kondo effecto'tl). 

The purpose of this paper is to discuss models of cou- 
pled impurities which can be tuned between regimes 
with different ground state spin. The spin quantiza- 
tion then requires that these regimes are separated by 
one or more boundary quantum phase transitions as the 
inter-impurity couplings are varied. Phases with a non- 
zero ground state spin are obtained in the situation of 
underscreeningQ, i.e., if a spin of size S interacts with K 
conduction electron channels, and 25' > K, leading to a 
residual free spin of size {S — K/2) as T ^ 0. A general 
two-impurity Kondo model necessarily showing a quan- 
tum phase transition then consists of two spins of size Si 
and 5*2, with a direct exchange coupling /, and a total of 
K screening channels, satisfying 2(S'i-f 5*2) > K. 

Related Kondo- and Anderson models have been dis- 
cussed recently in the context of multi-level or coupled 
quantum, dots with nearly degenerate singlet and triplet 
Ievelg3il3. Several authors have considered effective two- 
impurity models near the singlet-triplet degeneracy with 
two screening channelaij which arise in vertical multi- 
level dots; such a situation corresponds to the physics of 
the usual two-impurity Kondo model not showing a phase 
transition in the absence of particle-hole symmetryB. 
The case of a single screening channel and the associ- 
ated limiting situations of a singlet ground state and 
an undpscreened Kondo model have been mentioned in 
Refs.tjO, but the nature of the transition and the associ- 
ated quantum-critical behavior have not been examined 
to date. 

We note that a two-impurity jpodel with 2(S'i-|-S'2) < K 
can give rise to overscreeningQ in a channel-symmetric 
situation. Overscreening leads to a non-trivial interme- 
diate coupling fixed point already for a single impurity 
and is a very interesting topic on its own, which shall, 
however, not be discussed here. Remarkably, the two- 
impurity two-channel Kondo model (K = 4), has been 
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a) S = 1/2 b) S = 1 




FIG. 1. Kondo models discussed in the text. The 
solid/empty dots represent conduction electron orbitals and 
Kondo impurities, respectively. Ji,2 denote the Kondo cou- 
plings, I is an explicit inter-impurity exchange, i-aj) Two spins 
1/2, one channel, b) Two spins 1, two channelscj. 

solved exactly by Georges and Sengupta, andj-ahows a 
continuous line of non-Fermi liquid fixed pointstj. 

Two spins 1/2, one channel. The simplest realiza- 
tion of the above proposal is a "double Kondo impurity" 
(Fig. pi), i.e., a Kondo model of two spins 1/2 coupled to 
a single orbital of conduction electrons, with the Hamil- 
tonian H = i?band + -^int, 

ilint = JlSl -80 + J2S2 -So+ISi -82 (1) 

and iJband = J2ka ekCk^CkQ in standard notation, Sg = 
X)kk'a/3 '^laa'^ajscu'is IS the conduction band spin operator 
at the impurity site Tq = 0, and 81^2 denote the impurity 
spin operators. The Kondo couplings Ji, J2 are taken 
to be antiferromagnetic. For the numerical calculations, 
we will use a conduction band with a constant density of 
states, p{e) = po = 1/(2I>) for |e| < D = 1. 

The inter-impurity coupling / obviously drives a quan- 
tum phase transition: For large antiferromagnetic /, both 
impurities form a singlet with zero total spin, whereas 
large ferromagnetic / locks the impurities into an 5* = 1 
entity, which is underscreened by the single channel of 
conduction electrons, leading to a doublet ground state 
with residual spin 1/2. 

A second observation is that the symmetry w.r.t. to 
the two Kondo couplings plays a crucial role. Re- writing 
Hint (0) in terms of bond boson operators {s, 
with Si 2 = i^sH^ ± tjjS — ieai3-yt'jjtj)/2, one finds 

-ffint = I {tita - 3s'^s)/4 - ( Ji + J2) iea/S^SoiJji^ 

+ iJi-J2)s^{stl + h.c.). 

In the symmetric case, Ji = J2, there is no mixing be- 
tween impurity singlet and triplets, and the state formed 
by the inter-impurity singlet and a decoupled conduc- 
tion band is an exact eigenstate of H. The "rest" of 
the Hilbert space (triplet plus conduction band) forms 
an underscreened S — I Kondo problem. Therefore, for 
Ji ~ J2 the model shows a simple level crossing at a crit- 
ical Ic which marks a first-order transition between the 
singlet and the underscreened doublet state. 

Turning to the asymmetric case, Ji > J2, we start 
with results obtained using tiie Numerical Renormaliza- 
tion Group (NRG) techniqudlj. In Fig. ^ we display the 
temperature dependence of Tximp for fixed Ji_2 and dif- 
ferent values of /, where Ximp is the total impurity contri- 
bution to the uniform susceptibility. Let us first discuss 
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FIG. 2. NRG results0 for T^imp (its = 1) of the dou- 
ble-impurity model, a) Large asymmetry: Ji=0.1, J2 = 0.05. 
The thick curve corresponds to 7 = 0.00181527408 ~ Ic, the 
curves close to Ic are for 10^^(1 — Ic) = —3, -1-2, -1-3, +5, the 
others for |J - Ic\ = 10"", 10"*, 10"^ b) Small asymme- 
try: Ji=0.4, J2= 0.400001. Here h ~ 0.0788628547, curves 
arc displayed for 10^" (J-/c) = +0.5, +1, +1.2, +1.5, +7.5 and 

\i -ic\ = lo-^...,lo-^ 

the case of large asymmetry. Fig. |^a. We can identify 
Ic, with all curves for I > Ic leading to Tximp ^ 
for r — > (singlet), whereas all curves for I < Ic show 
Tximp 1/4, consistent with a residual spin 1/2. The 
topmost curves (/ ^ Ic) allow an identification of the 
Kondo temperature for the S = I Kondo problem, here 
we have Tk — 10~^^D (small values of Ji^2 have been 
chosen to illustrate the crossover behavior). The curves 
for I < Ic approach the fixed point very slowly, indi- 
cating an operator which is marginally irrelevant. For 
/ values near Ic and T 3> Tk, there is a plateau at 
Tximp — 1/2. For exchange couphngs / larger than but 
close to Ic, the system first flows to the fixed point with 
^^Ximp = 1/4, and eventually crosses over to the singlet 
ground state (for T/D < 10"^^ in Fig. ||a). The corre- 
sponding crossover scale T* is found to depend exponen- 
tially on the distance to the critical point, (I — Ic) - there 
is no power law which would characterize a second-order 
transition. Even for very small \I — Ic\ the NRG does 
not show any fixed point behavior which could be distin- 
guished from one of the stable fixed points, i.e., there is 
no sign of a third, infrared unstable, fixed point which 
could possibly correspond to the critical point at / = /c- 
We have analyzed the structure of the NRG fixed 
points, and found that for I > Ic the system consists of 
a singlet decoupled from the band ((Si • S2) = —3/4), 
whereas for all other values including Ic there is one 
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Kondo-screened spin 1 /2 and one residual decoupled spin 
1/2. Close to / = Ic, the NRG flow is identical to 
the one in the single-impurity Kondo model near the 
Kosterlitz-Thouless (KT) transition at J = 0. These 
results strongly suggest that the phase transition in the 
double-impurity model is of KT type, and the critical 
point is characterized by one decoupled spin- 1/2 degree 
of freedom. 

In the following, we give heuristic scaling arguments 
supporting this picture. Close to the critical point, we 
can assume that the effective inter-impurity coupling, 
Jeff = / -I- /rkky, is the smallest energy scale of the 
problem. In the spirit of poor man's scaling, both Kondo 
couplings will grow upon reducing the band cut-off, and 
the initial asymmetry between Ji and J2 increases under 
this process. The larger coupling, Ji, will become of or- 
der unity at an energy scale Tk ~ exp[— 13/(2 Ji)], in- 
dicating that the spin Si becomes screened. This screen- 
ing freezes the spin Si and the conduction electron spin 
So into a singlet, thus terminating the flow of J2 and 
effectively decoupling the spin S2 from the band. 

If we now switch on a small I^s <C Tr-, it will cou- 
ple the remaining spin via the screened singlet to the 
conduction band. This effective Kondo coupling can be 
estimated by second-order perturbation theory to be of 
order IcsD/Tk- Therefore, ferromagnetic Ics flows to 
zero as known from the ferromagnetic Kondo effect, and 
one ends up in a weak-coupling regime for the second 
spin, with a doublet ground state. On the other hand, 
antiferromagnetic I^s grows logaiithmically and leads to 
screening of the second spin 1/2LD. But, in contrast to 
a single-spin Kondo effect, this flow is cut-off when the 
effective coupling of the second spin becomes comparable 
to the singlet binding energy of the first spin Si, Tk, due 
to the indirect nature of this exchange. It is immediately 
clear what happens at this scale: The two spins Si and 
S2 acquire a large direct antiferromagnetic coupling, i.e., 
they form a singlet among themselves and effectively de- 
couple from the band. Therefore, this regime of weak 
antiferromagnetic /off is continuously connected to the 
one of strong antiferromagnetic / where the two spins 
form an inter- impurity singlet from the outset. We can 
identify I^s with {I — Ic), and the crossover scale T* near 
the critical point is given by T* ^ Dexp{—TK/[I — Ic])- 

In the case of maximal asymmetry, J2 = 0, the picture 
developed above is most transparent. Here, no RKKY 
interaction is generated, and the critical coupling value 
is obviously Ic — 0- For J2 > 0, we find Ic ~ J1J2/D > 
since the RKKY coupling is generically ferromagnetic be- 
cause the two antiferromagnetic Kondo couplings prefer 
a parallel alignment of Si and S2. Similady, J2 < (but 
Ji > 0) gives a KT transition at Ic < O.EI 

Summarizing, the critical point in the double-impurity 
problem is defined by the coupling value where one spin 
completely decouples from the band, and the other spin 
is Kondo-screened. The boundary quantum phase tran- 
sition is of the same type as in the single-impurity Kondo 
model when the bare coupling is tuned through zero, i.e., 
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FIG. 3. Schematicphase diagrams of the double-impurity 
Kondo model (Fig. hia). The ScS values in the different 
regimes are defined through the Curie term in the total impu- 
rity susceptibility, Tximp = Scff (S'cfi -I- 1) /3. a) For T — Q and 
fixed Ji. The solid line is a line of KT transitions, the open 
dot is the first-order transition for Ji = J2. b) For T > and 
fixed Ji,2 with not too small asymmetry (Ji — J2). The solid 
dot is the KT transition, the dashed lines denote crossovers. 

of the KT universality class0lii@. Close to the critical 
point for I > Ic, we find a new exponentially small energy 
scale T* within the Kondo regime - leading to two expo- 
nentially small scales, with T* <^ Tk - and the crossovers 
can be understood as a two-stage Kondo effect. Such a 
two-stagfi behavior was so far only known in two-channel 
modelsO'El; remarkably, here it is realized with a single 
conduction electron channel. 

We emphasize, however, that the stable fixed point for 
I > Ic is not given by two Kondo-screened impurities, 
but by an inter-impurity singlet. This has immediate 
consequences for the conduction electron phase shift. For 
ferromagnetic /cff {I < Ic), there is one screened spin 1/2 
leading to a phase shift of 7r/2, whereas for I > Ic the 
low-T limit of the phase shift is zero. 

It is also interesting to discuss the limit of small asym- 
metry of the Kondo couplings. As shown above, the 
singlet-doublet transition is of first order in the sym- 
metric case. Therefore, very small asymmetry introduces 
another low-energy scale, Tas oc | Ji — J2I, and the KT be- 
havior will be visible only for T ^ Tag. We have verified 
this behavior numerically, and sample results are shown 
in Fig. |b. Here, T^^-W'^^D, Tk^10~^D, and a new 
regime (compared to Fig. ||a) occurs for Tas < T < Tk 
where effectively level crossing between the singlet and 
the underscreened doublet is observed. Our numerical 
results are consistent with Tas ^ TkD\Ji — J2|/(Ji J2) ^ 
\Tki - Tk2\, with Tki,2 ~ Z?exp[-D/(2Ji,2)]. 

Schematic phase diagrams inferred from both the nu- 
merical results and the above arguments are shown in 
Fig. ||. AtT — 0, Fig. we have a line of KT transitions 
terminating at a first-order point for Ji = J2. The finite- 
T crossovers are escpecially interesting. Fig. ||b: Above 
Tk, the physics is dominated by the crossing of the S ~ 
and S — 1 impurity levels - in particular, near I — Ic 
the two spins 1/2 fluctuate independently, leading to the 
plateau at Tximp ^ 2 x S{S +l)/3 with 5" = 1/2. In con- 
trast, below Tk the behavior is determined by the KT 
quantum phase transition at Ic. 
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Two spins 1, two channels. We briefly want to com- 
ment on a different model, first considered in Ref. |2^, 
which consists of two spins 1 coupled to a total of two 
screening channels — i.e., two underscreened S = 1 
Kondo impurities, Fig. The ground state spin can be 
tuned between zero and one, and therefore there must be 
a quantum phase transition as the coupling / is varied. 
Ref. |2| did not consider such a transition- b.ut focussed 
exclusively on a possible phase transitiorijij associated 
with a critical non-Fermi liquid fixed point and a jump 
in the phase shift. It was found that such a transition 
does not exist for spins 1. 

In analogy to the discussion above, we expect that at 
sufficiently low temperatures and les 0, each impurity 
realizes an underscreened S = 1 Kondo effect, i.e., a spin 
1/2 remains unscreened at each of the impurity sites. 
Turning on a ferromagnetic /eft the two residual spins 
couple to an unscreened spin 1. In contrast, IcS > 
promotes an effective antiferromagnetic coupling, both 
among the two residual spins as well as to the conduction 
electrons. However, the energy gain due to singlet for- 
mation with the conduction band is exponentially small 
in Ics, whereas the energy of the inter-impurity singlet 
is linear in Jeff. Therefore, singlet formation between the 
two residual spins 1/2 is favorable, and this singlet will 
be effectively decoupled from the band. 

It follows that the boundary quantum phase transition 
in this model associated with the change in the ground 
state spin is always of first order, although it will show 
non-trivial finite-T crossover behavior due to residual fer- 
romagnetic Kondo couplings. Interestingly, this transi- 
tion does not affect the conduction electron phase shift, as 
the phase shift is determined by the two Kondo-screened 
spin- 1/2 degrees of freedom only, regardless of whether 
the two others form a decoupled singlet or triplet. 

Applications of the double-impurity model. Two impu- 
rities coupled to a single orbital may be hard to realize 
directly in metals, however, the double-impuritji model 
is obtained from the usual two-impurity modelcl in the 
limit of small inter-impurity distance R. In this limit, the 
screening channel associated with the odd linear combi- 
nation of conduction electrons is weak, and the physics of 
the double-impurity model can be observed over a wide 
temperature range. Experimentally, such a situation can 
be achieved with magnetic ad-atoms on a metal surface 
manipulated by scanning tunneling microscope (STM) 
techniquea23. 

Other applications can be-constructcd using coupled 
or multi-level quantum dotsE3. In the former case, each 
dot represents a single "impurity" spin, and the inter- 
impurity interaction / can be varied either by tuning 
the inter-dot tunneling rate or the direct exchange in- 
teraction. For multi-level dots, we point out that the 
single-channel situation considere d her e appears to be 
the generic case for lateral devicesE3E3 - in contrast to 
a two-channel model appropriate for vertical dots. In 
fact, the results of Ref. |2^ indicate that the KT transi- 
tion discussed for the double-impurity model can indeed 



be realized in a multilevel dot. 
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